Nonlinear excitations in CsNiF 3 in magnetic fields perpendicular to the easy plane 



CO 

o 
o 

CN 



o 

CO 

-4-* 



T3 

C 

O 

CN 
> 
O 

cn 

OS 

o 

cn 

o 

-I— > 



i 

c 

o 
o 



x 

S3 



M. Orendac, 1 A. Orendacova, 1 E. Cizmar, 2 J.-H. Park, 2 A. Feher, 1 S. J. Gamble, 2 S. Gabani, 3 K. Flachbart, 3 

J. Karadamoglou, 4 M. Poirier, 5 and M. W. Meisel 2 
1 Faculty of Science, P. J. Safdrik University, Park Angelinum 9, 041 54 Kosice, Slovakia 
2 Department of Physics and Center for Condensed Matter Sciences, University of Florida, Gainesville, FL 32611-8440, USA 
3 Institute of Experimental Physics, Slovak Academy of Sciences, Watsonova 47, 043 53 Kosice, Slovakia 
4 Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), EPFL, 1015 Lausanne, Switzerland 
5 Centre de Recherche en Physique du Solide, Departement de Physique, Univerite de Sherbrooke, J1K 2R1 Quebec, Canada 

(February 2, 2008) 

Experimental and numerical studies of the magnetic field dependence of the specific heat and 
magnetization of single crystals of CSN1F3 have been performed at 2.4 K, 2.9 K, and 4.2 K in 
magnetic fields up to 9 T oriented perpendicular to the easy plane. The experimental results 
confirm the presence of the theoretically predicted double peak structure in the specific heat arising 
from the formation of nonlinear spin modes. The demagnetizing effects are found to be negligible, 
and the overall agreement between the data and numerical predictions is better than reported for 
the case when the magnetic field was oriented in the easy plane. Demagnetizing effects might play a 
role in generating the difference observed between theory and experiment in previous work analyzing 
the excess specific heat using the sine-Gordon model. 
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I. INTRODUCTION 

Significant developments in solid state physics have oc- 
curred when nonlinear effects have been taken into con- 
sideration. Low dimensional magnets have proven to be 
excellent model systems for these studies since long-range 
ordering is suppressed, whereas short-range ordering, if 
sufficiently pronounced, gives rise to nonlinear behavior. 
In this respect, spin chains garner special importance, 
and the interest in these one-dimensional (ID) magnets 
was piqued by Mikeska's prediction of the contribution 
of thermally excited solitons to the dynamics of a planar 
spin chain. 1 These excitations are expected to appear 
in a symmetry breaking magnetic field oriented in the 
easy plane of the chain. Since the compound CsNiF 3 is 
a good realization of an S = 1 Heisenberg planar ferro- 
magnetic chain, it has become one of the most intensively 
studied quasi-lD systems. Its magnetic properties in a 
magnetic field oriented in the easy plane have been inves- 
tigated by different theoretical approaches 2,3 and various 
experimental techniques. More specifically, the presence 
of sine-Gordon solitons was assumed in the interpreta- 
tion of the central peak in the energy spectrum of slow 
neutrons, 4 of the temperature and field dependences of 
the spin-lattice relaxation time studied by NMR, 5 and of 
the excess magnetic specific heat. 6 In addition, the sup- 
presion of the energy gap in the excitation spectrum of 
CSN1F3 in a magnetic field oriented in the easy plane was 
analyzed using quantum sine-Gordon field theory. 7 

Although the sine- Gordon model has been widely ap- 
plied in the analysis of the data obtained on several 
Heisenberg chain systems, the existence of sine-Gordon 
solitons has not been proven unambiguously. For exam- 
ple, the temperature and magnetic field dependences of 
the intensity in the central peak observed by neutron 



scattering in CSN1F3 is better described by the spin- wave 
prediction in comparison with a model assuming scatter- 
ing from a dilute gas of solitons. 8 Furthermore, reason- 
able quantitative agreement between the experimentally 
measured peak of the excess specific heat and the cor- 
responding theoretical prediction could be obtained only 
after an artificial renormalization of the soliton energy, 6 
and this aspect has been the subject of considerable 
debate. 9 ' 10 It is noteworthy that a similar situation was 
found for the 5 = 5/2 Heisenberg planar ferromagnet 
TMMC. 11 In other work, optical studies were initially 
interpreted in terms of magnon-soliton scattering, 12 but 
were subsequently shown to be explainable by multi- 
magnon process. 13 Recent theoretical progress in calcu- 
lating the magnon dispersion and thermodynamic prop- 
erties of CSN1F3 using a transfer matrix renormalization 
group algorithm 3 has allowed for the comparison of the 
experimental data with the exact theoretical predictions. 
Nevertheless, using generally accepted values of the ex- 
change coupling J, the single- ion anisotropy D, and the g 
factor for CsNiF3, only qualitative agreement was found 
for the magnetic field dependence of the specific heat. 

Since solitons are predicted to appear when the mag- 
netic field is oriented in the easy plane, less atten- 
tion has been devoted to the situation when the field 
is oriented perpendicular to the easy plane. Apart 
from the experimental study of the temperature depen- 
dence of the susceptibility 14 and the field dependence of 
the magnetization 15 for such a field orientation, other 
work focused predominantly on mapping the magnetic 
phase diagram obtained from the anomalous behavior of 
the elastic constants. 16 In addition, a theoretical treat- 
ment addressed the field induced instability of the three- 
dimensional (3D) antiferromagnetically ordered phase at 
T < Tjv = 2.7 K. 17 
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In fact, the situation is significantly richer than ini- 
tially expected due to the formation of nonlinear spin 
modes when the magnetic field is orientated perpendic- 
ularly to the easy plane. 3 More specifically nonlinear ef- 
fects lead to the formation of a characteristic double peak 
in the magnetic field dependence of the specific heat in 
the vicinity of a critical field B c at which the ground 
state is changed. The existence of the double peak has 
already been predicted in S = 1/2 XY chains as well 
as S = 1/2 Heisenberg chains with Ising anisotropy. 18 ' 19 
The proposed behavior of the specific heat significantly 
differs from that obtained in the dilute-magnon approx- 
imation, where a single broad maximum at B c is pre- 
dicted. The motivation of our experimental study of 
CsNiF 3 was to search for the existence of the predicted 
double peak structure in the field dependence of the spe- 
cific heat and to examine the thermodynamic response 
at several temperatures near Tjv. As we will present in 
the following sections, our results confirm the existence of 
the theoretically predicted double peak structure. More- 
over, the numerical analysis of our results indicates that 
the previously reported discrepancies between the excess 
specific heat data and the corresponding theoretical pre- 
dictions may be partially attributed to demagnetization 
effects. 



II. EXPERIMENTAL DETAILS 

The single crystals of CSMF3 used for our work were 
light green and transparent throughout their volume, al- 
beit with some internal cleavage planes noticeable. Due 
to the highly hydroscopic nature of the material, all 
samples were kept in paraffin oil until they were ready 
to be mounted. For the specific heat study, a single 
crystal of mass 8.44 mg and approximate dimensions 
7 x 0.5 x 0.5 mm 3 was used. This specimen was attached 
to a simple homemade microcalorimetcr using Epotck 
E4110 silver epoxy. The calorimeter consisted of a sap- 
phire disk that supported a heater and a thermometer, 
both made from RuC>2 chip resistors. The thermometer 
was calibrated against a commercial Cernox thermome- 
ter (Lakeshore CX-1030CD) which was located near the 
calorimeter. The position of the calorimeter was fixed by 
a holder made of Vespel SP-1 that was needed to pre- 
vent motion in the magnetic field. The holder and the 
manganin wires from the thermometer and heater served 
as the thermal link to the stable reservoir. At 4.2 K, the 
thermal conductivity of the manganin wires (0.15 /iW/K) 
represented only about 5% of that of the vespel holder 
(3 /iW/K). Consequently, the resultant value of the total 
thermal conductivity can be considered to be field inde- 
pendent. The resistance of the thermometer was mea- 
sured using a Wheastone bridge, with a PAR 124A lock- 
in amplifier serving as a null detector. While monitoring 
the cooling curve of a measurement, the analog output 
of the lock-in was read by a HP 3457A multimeter. The 



same device was used for measuring the voltage on the 
heater during the heating step, and the current was pro- 
vided by a Keithley Model 220 current source. Conse- 
quently, the temperature and magnetic field dependences 
of the heater resistance did not influence the accuracy of 
the measurements. In the measuring cycle, the relative 
change of the temperature of the sample ranged from 3% 
to 5%. The specific heat was calculated from the ratio of 
the characteristic cooling relaxation time and the ther- 
mal conductivity of the link. We estimate the overall 
accuracy of the measurement to be better than 5%. 

For the estimation of the correction due to demagne- 
tizing effects, the field dependence of the magnetization 
was studied up to 7 T at several selected temperatures 
using a Quantum Design SQUID magnetometer. For 
these measurements, a single crystal, with dimensions 
similar to those used for specific heat studies, was placed 
in a polyethylene vial that was held by a straw. The 
background contribution from the vial and straw was in- 
dependently measured to be less than 1% of the total 
signal and was subtracted. 

III. RESULTS 

For our first step, the specific heat of CSN1F3 was stud- 
ied from 2 K to 5 K in zero magnetic field. These re- 
sults were compared with the data reported by other 
workers, 20 and the comparison is presented in Fig. 1. 
Both sets of data are characterized by a A-like anomaly at 
about 2.7 K, indicating the presence of long-range order- 
ing. However, the larger peak observed in our experiment 
suggests that our single crystal specimen was of higher 
quality than the one used in the other work. 20 The excel- 
lent agreement of the critical temperatures confirms that 
the silver epoxy did not deteriorate the bulk properties 
of our sample. The difference between both data sets can 
be described by the equation 

C(T) = aT + bT 3 , (1) 

when a = 0.9 ± 0.3 J/K 2 and b = 0.88 ± 0.03 J/K 4 . This 
difference is attributed to the specific heat of the addenda 
created predominantly by the vespel holder and the sap- 
phire substrate. In addition, since the possible amount 
of paramagnetic impurities in the silver epoxy is at most 
0.001%, 21 the specific heat of the addenda should have a 
negligible magnetic field dependence. The fact that the 
lattice contribution of CsNiF 3 was not subtracted from 
the total heat capacity does not compromise the evalua- 
tion of the magnetic field dependence of the excess spe- 
cific heat, the quantity of primary interest in this study. 

The magnetic field dependence of the specific heat was 
studied at 2.4 K, 2.9 K, and 4.2 K, in fields up to 9 T. 
As mentioned previously, the weak magnetoresistance of 
the heater is simple to accommodate and has no direct 
influence in a thermal relaxation study. Similarly, the 
weak magnetoresistance of the thermometer is negligible 
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since it docs not influence the thermal relaxation time 
of the calorimeter. On the other hand, the data were 
corrected using the known magnetoresistance of the Cer- 
nox thermometer. 22 As expected, the correction was most 
significant at T = 2.4 K in magnetic fields higher than 
6 T, where it represented 3.5% of the measured heat ca- 
pacity value. The resultant magnetic field dependence of 
the total specific heat is presented in Fig. 2. The A-like 
anomaly observed at 2.4 K in a field of about 2 T corre- 
sponds to field-induced long-range ordering. The values 
of the critical temperature and field agree well with the 
magnetic phase diagram obtained from ultrasonic inves- 
tigations of specimens from the same batch of crystals. 16 



IV. DISCUSSION 

For the investigation of the effect of the nonlinear ex- 
citations on the equilibrium thermodynamic properties, 
most of the attention has concentrated on the excess 
isothermal specific heat 



AC(T,B) = C(T,B) - C(T,0) , 



(2) 



studied as a function of magnetic field B at a temperature 
T. The numerical predictions 3 of the excess specific heat 
for CsNiF 3 in a magnetic field oriented perpendicular to 
the easy plane were derived from a ID model described 
by the S = 1 Hamiltonian 

n = -jJ2 s i- s i+i + dJ2( s * z ) 2 + g^ B Bj2s t z , 

i i i 

(3) 

using the values J/ks = 23.6 K, D/ks = 8.25 K, and 
g = 2.13. The numerical simulations were performed for 
temperatures 2.4 K and 4.2 K, the former being below 
the critical temperature of the real material. The results 
of the numerical simulations suggest that the predicted 
double peak structure in the excess specific heat becomes 
more pronounced at lower temperatures. Consequently, 
in the corresponding experimental study, the influence of 
interchain correlations should be considered in the criti- 
cal region. 

The ratio of intrachain and interchain exchange cou- 
pling constants (J/ J' w 500) makes CSN1F3 one of the 
best representatives of a quantum spin chain. 23 Indeed, 
the long-range ordering observed at 2.66 ± 0.01 K in this 
study is manifested as a small spike located on the broad 
maximum appearing due to the short-range correlations. 
In addition, the amount of entropy removed by the spike 
itself (w 0.05 J/(K mol)) is by far less than 1% of the 
total amount of entropy for the S = 1 system, namely 
9.13 J/(K mol). In such a situation, it appears reason- 
able to approximate the contribution of the interchain 
coupling to the total specific heat by fitting the temper- 
ature dependence of the specific heat using data below 
and above the spike, i.e. outside the critical region. The 



result of this approximation is shown as the solid line 
in Fig. 1, and the enhancement of the specific heat in 
the critical region due to the interchain correlations is 
apparent (see inset in Fig. 1). Such an enhancement 
generates a deviation from the corresponding theoretical 
predictions, 3 which were made for a purely ID system. 
Consequently, when comparing the experimental excess 
specific heat data obtained at 2.4 K with the theory, the 
specific heat value calculated from the aforementioned 
fit at 2.4 K was taken as the reference value C(T, 0) in 
Eq. (2). As a result, if the experimental specific heat data 
are represented by C(T,B) in Eq. (2), then the quan- 
tity AC should tend to a nonzero value in the zero field 
limit. In other words, in the critical region, the inter- 
chain coupling is responsible for the finite value of the 
excess specific heat. The comparisons of the experimen- 
tal excess specific heat data with the numerical predic- 
tions are presented in Fig. 3. As expected at 2.4 K, the 
agreement between the theory and experiment is greatly 
improved in magnetic fields greater than « 3 T, which is 
high enough to overcome the influence of the long-range 
correlations and decompose the system into independent 
chains. This difference between the experimental data 
and theoretical predictions persists, but is smaller, at 
2.9 K, and it disappears at 4.2 K. It should be noted that 
critical magnetic field observed at 2.4 K is consistent with 
the phase diagram obtained from the ultrasonic study 
which used specimens from the same batch of crystals. 16 
Furthermore, the experimental data reproduce the pre- 
dicted positions of the minimum values of AC, as well as 
the shift of the position of the second maximum towards 
lower fields with decreasing temperature. Although some 
quantitative differences persist, it should be stressed that 
the overall agreement between the numerical predictions 
and experimental data is better than that achieved when 
the magnetic field is oriented in the easy plane. 3 

Several potential reasons might explain the observed 
differences. For example, quantitative comparisons be- 
tween the experimental data and the theoretical predic- 
tions should take demagnetizing effects into account. The 
correction of the specific heat due to demagnetizing ef- 
fects can be calculated as follows. 24 The internal mag- 
netic field is written as 



Bi = B„ - NM 



(4) 



where Bi and B e stand for the internal and external 
fields, while M and N denote the magnetization and de- 
magnetizing factor. The specific heat of the sample may 
then be written as 



Cb^B^^CbA^B,) 



[NT(dM/dT) 2 Be ] 
[1 - N(dM/dB e ) T ] 



(5) 



The needle shape of the sample, which was oriented with 
its long axis parallel with the applied external magnetic 
field, should lead to a small value of the demagnetizing 
factor. Indeed, for the dimensions of our single crystal 
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specimen, the average demagnetizing factor is approxi- 
mately 0.04 in SI units. In order to evaluate (dM/dB e )T, 
the field dependence of magnetization was investigated 
up to 7 T at 2.4 K, 2.9 K, and 4.2 K. For clarity, only 
the results obtained at 2.4 K and 4.2 K are presented in 
Fig. 4. For all temperatures, the magnetization is charac- 
terized by an increase up to the critical field, followed by 
the tendency to saturate at higher fields. It should also 
be noted that, for the magnetization data corrected for 
demagnetizing effects, some deviation from the numerical 
predictions appears. The values of the term {dM/dT) Bc 
were evaluated with the help of the additional field de- 
pendences of the magnetization studied at temperatures 
about 2% higher and lower than the aforementioned ones. 
Then (dM/dT) Bc for 2.4 K, 2.9 K, and 4.2 K were ap- 
proximated by the corresponding differences of the mag- 
netization at the neighboring temperatures. The results 
of this analysis are presented in Fig. 5. The remarkable 
feature of the (dM / 8T)b c quantity is its nontrivial field 
dependence in the region of low fields, where it is most 
pronounced at 2.4 K and also observed at 2.9 K, but 
is completely absent at 4.2 K. Taking into account that 
both 2.4 K and 2.9 K belong to the critical region of 
temperatures for CSN1F3, the observed behavior can be 
attributed to the effect of interchain coupling. It should 
be noted that, for 2.4 K and 2.9 K, the value of the field, 
where the anomalous behavior of (dM/ dT) b c disappears 
(«3 T), agrees very well with that for which AC starts 
to follow the theoretical prediction for a pure ID system. 
Such agreement supports the suggestion that « 3 T is 
sufficient to suppress the 3D effects near Tn. After eval- 
uating (dM/dB)x and (dM/dT)s c , the correction to the 
specific heat due to the demagnetizing effects was calcu- 
lated. However, it turns out that the correction does not 
exceed 0.1% of the uncorrected specific heat value, thus 
demagnetizing effects can not be responsible for the ob- 
served differences. This fact contrasts with the situation 
when magnetic field is oriented in the easy plane, where 
such a correction was proven to be significant. 25 

Naturally, tuning the parameters J, D, and g might 
improve the agreement between the numerical predic- 
tions and the experimental results. Consequently, the 
susceptibility, 14 magnetization, 15 and present specific 
heat data were reanalyzed using the transfer matrix 
renormalization group technique. Since no evidence for 
the change of exchange coupling constant was found in 
the previous analysis, 3 the value J/k B — 23.6 K was 
adopted in the recalculation. In the new set of param- 
eters, the 'standard' value of D/ks = 9 K was cho- 
sen, whereas g = 2.19 was obtained from fitting the 
susceptibility 14 and magnetization 15 with the new D 
value. It should be noted that such a choice of new 
D and g parameters does not significantly change the 
magnitude of the critical field, B c = D/gfiB ~ 6 T as 
suggested by the specific heat data. However, as can be 
seen in Fig. 3, the new predictions are shifted towards 
larger fields and the values of AC are increased, thereby 
making the agreement with the experimental data worse 



than for the original set of parameters. Further increas- 
ing of the D value leads to unsatisfactory fitting of the 
susceptibility and magnetization. Consequently, tuning 
the parameters in the framework of a pure ID model does 
not suppress the differences between the numerical pre- 
diction and the experimentally measured excess specific 
heat. 

Finally, the possible effect of the tilting of the c-axis 
with respect to the direction of the magnetic field may 
be considered. It is not straightforward to quantitatively 
evaluate this effect, but it can be expected to influence 
the experimental results. For example, the susceptibility 
below 8 K is about five times higher for magnetic field 
oriented in the easy plane than that for the perpendicular 
orientation. 3 Although the magnitudes of AC are compa- 
rable for the different field orientations, the correspond- 
ing field dependences are completely different. Thus, tilt- 
ing may nontrivially contribute to the differences between 
the numerical predictions and experimental specific heat 
data. 



V. CONCLUSIONS 

The experimental study of the magnetic field depen- 
dence of excess specific heat performed at 2.4 K, 2.9 K, 
and 4.2 K has confirmed the theoretically predicted dou- 
ble peak structure reflecting the nonlinear behavior of the 
system. The excess specific heat data reproduce the pre- 
dicted shift of the second maximum towards lower fields 
with decreasing temperature. It was found that, simi- 
larly as for the orientation of the magnetic field in the 
easy plane, it is not possible to fit all available thermo- 
dynamic data using a pure ID model with a single set of 
parameters. Although tilting of the sample may be an 
alternative explanation of the persisting deviations, de- 
magnetizing effects arc determined to be negligible. This 
conclusion contrasts with the one made in Ref. 6, where 
the magnetic field was oriented in the easy plane and 
where the correction for the shape dependence had a pro- 
nounced effect. However, even with the incorporation of 
the demagnetization correction, the agreement between 
the excess specific heat data 6 and the corresponding ex- 
act numerical predictions 3 is worse than in our case. In 
Ref. 6, the estimation of (dM/dT) B<1 from the field de- 
pendence of the magnetization 26 might not be sufficiently 
accurate due to the restricted amount of data that is 
available. Indeed, in situations where the demagnetizing 
effects are pronounced, more detailed mappings of the 
field and temperature dependences of the magnetization 
were performed. 27 Naturally, a question may arise about 
the influence of demagnetizing effects in the specific heat 
studies of other soliton - bearing systems. A considerable 
amount of theoretical effort 28-30 has been developed to 
understand the differences and to improve the agreement 
with the specific heat data. 11 ' 31 However, to the best of 
our knowledge, the data themselves were not corrected 
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for demagnetizing effects. Consequently, a future exper- 
iment might focus on clarifying the role of the demag- 
netizing effects in CSN1F3 as a soliton - bearing system. 
Furthermore, since the predicted double peak structure 
seems to be a generic feature of the nonlinear behavior 
in classical and quantum spin chains, future experimen- 
tal effort may focus on measuring AC in other related 
materials. 
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FIG. 1. The temperature dependence of total heat capacity 
of CsNiF3 together with the addenda (squares) is compared 
to the data of magnetic heat capacity reported in Ref. 20 (cir- 
cles), after renormalization to the mass of our sample. The 
estimated contribution of the addenda for the present study 
is shown by the triangles. The heat capacity of our sample, 
i.e. the difference between the squares and triangles, is rep- 
resented by the crosses. A least square fit of the lattice con- 
tribution and the addenda is denoted by the dotted line, see 
Eq. (1). The solid line is obtained by fitting the data outside 
the critical region, see the text for a more detailed discussion. 
The inset provides an expanded view near the critical region. 
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FIG. 2. The magnetic field dependence of the total specific 
heat of CsNiF3 is shown at 2.4 K (triangles with lines as 
guides for the eyes), 2.9 K (squares) and 4.2 K (circles). 
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FIG. 3. The magnetic field dependence of AC, Eq. (2), 
at (a) T = 2.4 K, (b) T = 2.9 K, and (c) T = 4.2 K, 
are compared with theoretical expectations. The solid lines 
represent the numerical predictions when J/ks = 23.6 K, 
D/ks = 8.25 K, and g = 2.13, while the dotted lines are for 
J jk B = 23.6 K, D/k B = 9 K, and g = 2.19. 
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FIG. 4. Magnetic field dependence of magnetization of 
CsNiF 3 studied at 2.4 K (squares) and 4.2 K (circles). The 
solid and dotted lines represent the numerical predictions 
at 2.4 K and 4.2 K, respectively, when J/fcs = 23.6 K, 
D/k B = 8.25 K, and g = 2.13. 
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FIG. 5. Magnetic field dependence of (9M/9T)s c calcu- 
lated at 2.4 K (solid line), 2.9 K (dashed line), and 4.2 K 
(dotted line). 
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